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•Simple  analytic  approximate  expressions  have  been  derived  for  ! 
propagation  constants  in  the  earth-ionosphere  waveguide.  The 
expressions  are  an  extension  of  previous  work  to  a  more  genera, 
class  of  ionospheric-conductivi ty  height,  profiles.  Tin  oxpros; 
derived  here  allow  the  rapid  computation  of  ELF  phase  speeds, 
attenuation  rates,  and  excitation  factors  for  a  wide  r  ingo  of  i 
osphere  conditions  without  the  necessity  of  lengthy  lull-wave 
computer  calculations.  The  results  can  bo  annliod  to  Mo  tapi, 
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I .  INTRODUCTION 

In  two  recent  papers  (Greifinger  and  Greifinger  [1978] 
and  Greifinger  and  Greifinger  [1979],  hereinafter  referred  to 
as  Paper  I  and  Paper  II) ,  the  authors  derived  simple  analytic 
approximate  expressions  for  TEM  eigenvalues  for  ELF  propaga¬ 
tion  in  the  earth-ionosphere  waveguide,  assuming  that  the 
ionospheric  conductivity  profile  could  be  approximated  in  a 
certain  convenient  way.  These  expressions  make  possible  a 
rapid  calculation  of  ELF  phase  velocities  and  attenuation 
rates  without  the  necessity  of  lengthy  full-wave  numerical 
calculations.  Comparisons  were  made,  for  a  number  of 
ionospheric-conductivity  profiles,  between  approximate  and 
fall-wave  eigenvalues,  and  the  results  were  found  to  be  in 
excellent  agreement  in  all  cases.  More  recently,  Bannister 
[1979]  has  used  our  approximate  theory  and  the*  Wait  [1964, 
1970]  VLF  exponential  ionospheric-conductivity  profile  to 
determine  propagation  constants  for  ELF  daytime  propagation 
in  the  earth- ionosphere  waveguide.  He  showed  that  the  result¬ 
ing  values  of  ELF  attenuation  rate  and  phase  velocity  are  in 
excellent  agreement  with  measured  data  from  various  sources  in 
t.ie  frequency  range  from  5  to  200U  Hz. 

The  work  of  Bannister  clearly  demonstrates  the  validity 
and  calculational  utility  of  the  approximate  theory.  It  also 
adds  credence  to  the  belief  that  the  theory  should  be  equally 
applicable  to  PCA-disturbed  ionospheres,  for  which  comparable 
da* a  ao  not  yet  exist,  and  to  njclcar-dis turbed  ionospheres. 

As  such,  it  is  potentially  very  useful  for  making  rapid  and 
inexpensive  sensitivity  studies  of  ionospheric  propagation 
under  disturbed  conditions.  It  therefore  seems  useful  to 
extend  the  approximate  theory  to  a  more  ueneral  class  of 
ionospheric-conductivity  profiles  than  those  considered  in 
Papers  I  and  II.  The  conductivity  profiles  to  which  the 
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results  in  Papers  I  and  II  apply  will  appear  as  limiting  cases 
of  the  more  general  profile.  In  addition/  approximate 
analytic  expressions  will  be  derived,  within  the  framework  of 
the  theory,  for  the  excitation  factor,  which  is  an  important 
additional  measurable  quantity  in  ELF  propagation  experiments. 
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SUMMARY  OF  PREVIOUS  WORK 


II . 


In  Paper  I,  the  authors  developed  the  theory  for  an 
ionosphere  that  was  sufficiently  disturbed  that  there  was  no 
significant  penetration  of  the  electromagnetic  field  to 
altitudes  where  anisotropy  due  to  the  Earth's  magnetic  field 
had  to  be  taken  into  account.  It  was  shown  that  the  propaga¬ 
tion  constant,  under  these  conditions,  was  determined  by  four 
parameters  —  two  frequency-dependent  altitudes  and  the  local 
conductivity  scale  height  at  each  of  those  altitudes.  The 
lower  altitude,  denoted  by  h ,  is  the  altitude  at  which  the 
conduction  current  becomes  equal  to  the  displacement  current. 
The  higher  altitude,  denoted  by  h^,  is  where  the  reciprocal 
of  the  local  wave  number  becomes  equal  to  the  local  scale 
height  of  the  refractive  index.  If  the  conductivity  profile 
in  the  neighborhood  of  h  and  h.  can  be  adequately  annroxim.at 
by  a  local  exponential,  the  expression  for  the  TLM  o iqenva i uc 
is 
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where  ' g  and 
respectively 
velocity  is 


arc  the  local  scale  heights  at  hg  and  h^, 

From  the  definition  of  S  ,  the  relative  i huso 
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and  the  horizontal  attenuation  rate  in  decibels  per  mecjametor 


a  =  8.68  x  10 ^  -  Im (S  ) 
c  o 
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where  f  is  the  frequency  in  Hertz. 

In  Paper  II,  the  theory  was  extended  to  include  the  effects 
of  anisotropy.  For  the  sake  of  mathematical  simplicity,  it 
was  assumed  that  the  Earth's  field  was  vertical,  but  the 
results  actually  apply  over  a  rather  large  geographic  range. 
When  anisotropy  is  included,  the  parameters  which  determine 
tiie  propagation  constant  differ  for  daytime  and  nighttime 
ionospheric  conditions.  For  both  sets  of  conditions,  two 
parameters  which  enter  are  the  frequency-dependent  altitude 
Hq  at  which  tiie  conduction  current  parallel  to  the  magnetic 
fielu  becomes  equal  to  the  displacement  current,  and  the  local 
scale  height  of  the  parallel  conductivity  a  .  Under  daytime 
conditions,  two  additional  pairs  of  frequency-dependent 
altitudes  and  associated  scale  heights  appear  as  parameters. 

One  altitude  is  that  at  which  the  reciprocal  of  the  local  wan 
number  for  vertically  propagating  O  waves  becomes  equal  to 
the  local  scale  height  of  the  refractive  index  for  such  waves. 
The  associated  scale  height  is  the  scale  height  of  this  refrac¬ 
tive  index.  The  other  pair  of  parameters  are  the  corresponding 
quantities  for  vertically  propagating  X  (whistler)  waves.  If 
these  altitudes  arc  attained  in  a  region  of  the  ionosphere 
wave  jujj  |Cp|,  where  o  is  the  Hall  conductivity  and  ^  the 
E’edcrsen  conductivity,  the  two  pairs  of  parameters  become  iden¬ 
tical.  This  is,  in  fact,  the  case  over  a  substantial  altitude 
range.  Under  these  conditions,  the  approximate  expression  i or 
the  TEM  eigenvalue  is 


where  the  altitude  h^  and  scale  height  t  are  referred  to  the 
Hall  conductivity.  This  expression  is  slightly  different  from 
Eq.  (1),  its  counterpart  for  the  isotropic  ionosphere.  How¬ 
ever,  since  '0/hQ  <<  1  and  c1/h-L<<l,  the  value  of  SQ  is  the 
same  to  lowest  order  in  these  quantities. 

For  typical  nighttime  conditions,  a  sharp  reflecting 
E-region  bottom  may  be  encountered  before  the  reciprocal  of 
the  local  wave  number  becomes  equal  to  the  local  scale  height. 
Under  such  circumstances,  the  altitude  of  the  E-region  bottom 
replaces  h^  as  a  parameter,  and  the  local  wavelength  on  the 
E-region  side  of  the  bottom  replaces  the  scale  height  as  a 
parameter.  Under  these  conditions,  the  approximate  TEM 
eigenvalue  is  given  by 


where  =  e/c  and  n  is  the  index  of  refraction  for  X  waves 
just  inside  the  E-region. 

The  approximate  daytime  expressions  were  derived  under  (he 
assumption  that  the  conductivity  profile  in  the  neighboi  h< kks 
of  h^  could  be  approximated  by  a  local  exponential  over  an 
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altitude  range  of  a  scale  height  or  so.  The  approximate 
night  rime  expression  was  derived  under  the  assumption  of  a 
sharp  transition  in  the  ionospheric  conductivity  profile  he J ow 
which  the  local  wavelength  is  large  compared  with  the  local 
scale  height  and  above  which  it  is  small  compared  with  the 
local  scale  height.  There  is  an  important  class  of  ionospheric- 
conductivity  profiles  which  bridge  the  gap-  between  the  assumed 
daytime  conditions  and  nighttime  conditions.  These  are  condi¬ 
tions  where  the  altitude  h^  is  attained  in  a  region  of  rapid, 
but  not  discontinuous,  increase  of  conductivity  with  altitude . 
Under  such  circumstances,  the  conductivity  ro file  in  the  nei- 
borhood  of  the  altitude  h^  cannot  ho  accurate iy  r« presentee  by 
a  local  exponential.  In  tiro  following  se  .tier; ,  the  ay;  ro.\  im.ate 
theory  of  Papers  I  and  II  will  be  extended  to  deal  with  sues  a 
s i tua t ion . 


IT  I.  GUNLRALIZUD  CONDUCTIVITY  PROl'ILK 


In  this  section,  we  consider  a  more  General  ionospheric- 
conductivity  profile-  in  the  neiqhborhood  of  h.  than  those 
treated  in  Papers  I  and  II.  As  in  the  previous  work,  is 
defined  as  the  altitude  at  which 
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where  ^  is  the  conductivity  at  the  altitude  h^  am. 


1  do  1-1 


is  the  local  scale  height  at  this  altitude.  A  more  general 
conductivity  profile  in  the  vicinity  of  hj  can  be  written  as 


( z-h. ) / 

( 1  +  s )  c 

1  (z-h,)/' 

1  -  sc 


where  s  is  an  additional  parameter.  The-  scale  height 
related  to  the  other  parameters  through 


2  i  1  +  s) 


which  follows  from  let  (H)  air.  (9)  . 
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=  e 


(hl  -  hE}/  2 


whore  Eq.  (11)  defines  the  altitude  h£ .  The  profile  given  by 
Lq.  (9)  then  takes  the  form 


( z-h  ) / : 

(1  +  s)  e _ 

1  (z-h  )/' 

1  +  e  b  ^ 


(12) 


This  profile  is  approximately  exponential  up  to  altitudes  a 
scale  height  or  so  below  h  and  becomes  essentially  constant 

il< 

a  scale  height  or  so  above  h^ .  This  is  the  kind  of  transition 
which  typically  occurs  between  the  D-region  and  the  E-region, 
and  h  may  thus  be  thought  of  as  the  lower  boundary  of  the 
E-region.  If  h  >>h^,  the  profile  in  the  neighborhood  of 
is  a  simple  exponential  with  scale  height  • - •  ,  which  was 

the  daytime  profile  treated  in  Papers  I  and  II.  If  h  h^, 
the  E-rcgion  is  encountered  before  the  altitude  is  attained. 
If,  in  addition,  the  scale  height  •  is  sufficiently  small,  the 
transition  between  the  D-  and  E-regions  is  very  sharp,  which 
was  the  nighttime  profile  treated  in  Paper  II.  The  generalized 
conductivity  profile  given  by  Eq.  (9)  thus  contains  the  pre¬ 
vious  profiles  as  limiting  cases. 
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IV. 


DETERMINATION  OF  APPROXIMATE  TEM  EIGENVALUES 


As  in  Paper  II,  we  assume  that  the  Earth's  magnetic  field 
is  vertical.  In  the  QL  approximation,  propagation  in  the 
vicinity  of  h^  is  governed  by  two  uncoupled  wave  equations. 

In  terms  of  quantities  a  defined  in  Paper  II,  these  can  be 
written  as 


jZ 


2,2  , 

n,V 
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(13) 


where 


»2  *  -Jr  «*••>  •  W) 

The  upper  sign  corresponds  to  vertical  0  wave  propagation  and 
the  lower  sign  to  vertical  X  wave  (whistler)  propagation. 

Under  normal  daytime  conditions  at  the  upper  end  of  the  ELF 
band  and  for  sufficiently  disturbed  daytime  or  nighttime 
conditions  at  all  ELF  frequencies,  the  altitude  h^  is  attained 
in  a  region  where  the  ionosphere  is  nearly  isotropic,  so  that 


O  ' 

at  low  frequencies  and  normal 
frequencies,  h^  occurs  where 


H  *  °'  V'  -'0  and 


2  2 
n+  =  n_ 


Under  normal  daytime  conditions 
nighttime  conditions  at  all  ELF 
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so  that 


1 1 6 ) 
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JL!  _ 
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It  will  be  assumed  that,  whichever  case  applies,  the  relevant 
conductivity  profile  (oQ  or  .•■  )  is  given  by  Eq .  (9).  Both 

cases  may  be  treated  within  the  same  framework  by  writing 
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(1  +  s) 
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(z-hf) /’ 2 
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where 
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(18) 


for  the  isotropic  case  and 
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(19) 


for  the  anisotropic  case. 

The  wave  equation  may  be  solved  by  introducing  the  variable 

( z-hj^)  A. 2 

y  =  -se 

which,  with  Eqs.  (7)  and  (10),  transforms  Eq .  (13)  to 

..2  i". 

+  L  __1‘  +  _ I _ 

2  y  ty  4s  ( 1  +  s)  y  (y  -  1 ) 


-  0. 


(,'v) 


This  can  bo  recognized  as  the  generalized  hypergeometric 
differential  equation.  The  solutions  which  obey  the  radia¬ 
tion  condition  at  large  altitude  (y  v  - • )  are  (Magnus, 
Obcrhettinger ,  and  Soni  [1966]) 


1  +  2\-  (  ; 


C- 


where  F  is  the  generalized  hypergeometric  function,  and 


i  ( -  7 r)  /  2 
e _ 

2[s(l  +  s)P 

These  solutions  are  valid  down  to  an  altitude  a  few  scale 
heights  below  h^  or  hE ,  whichever  is  smaller.  At  these 
altitudes  (y  -  0)  ,  the  solutions  take  the  form  (Magnus, 
Obcrhettinger,  and  Soni  [1966]) 


*  C  ,  [  z  -  G  ( i-  ,  )  ] 


where  C,  are  constants,  and 


0  (  ;■  ) 


2'.  (1  +  ■  )  -  2,(1) 


-  +  log  s 


where  s  is  the  logarithmic  derivative  of  the  gamma  function 


The  eigenvalue  is 


s<  >  i  u  t  i  <  -nr. 
Wound  and 


winch  obey 
are  valid 
sol ut ions 


determined  by  matching  Kg.  ( )  to 
the  proper  boundary  conditions  it  tin 
up  to  altitudes  a  few  scale  heights  ai 
were  shown  in  Paper  1  I  to  have  tin  let 


i’aesi 


=  A  1(1  i  C)  z  -  S' 


where  A  and  C  are  constants.  For  z  >>  h^,  the  integrand  in 
Eq.  (25)  becomes  very  small,  and  little  error  is  made  by 
extending  the  range  of  integration  to  infinity.  Matching 
Eqs,  (23)  and  (25)  in  the  region  where  both  are  valid  then 
leads  to  the  condition 


2G  ( |3+)  G(E_) 

g  ( )  +  G(e_r 


which  is  the  generalized  expression  for  the  TF.fl  eigenvalue. 


V. 


COMPARISON  WITH  PREVIOUS  RESULTS 


The  approximate  expressions  derived  in  Papers  I  and  II  ear. 
be  recovered  as  limiting  cases  of  Eq.  (20.  In  those  papers, 
it  was  assumed  that  the  conductivity  profile  in  the  neigh¬ 
borhood  of  h^  coula  be  represented  by  a  local  exponential , 

i.  e.  , 
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iz~h0}/'Q 
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(2V) 


With  this  approximation,  the  integral  in  Eq .  (2f)  can  bo 

evaluated  analytically ,  giving 
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The  eigenvalue  can  thus  be  written  as 


f  ( 


+  ' ' 


h0  -  2 
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where 


20  (  4  )  f,(  _) 

f  ( '  +  '  •’  -  }  (TU-  JTgTj  ' 


The  limiting  forms  of  f  )  will  now  In  •  net  <  t  i  .  i . 

a.  Daytime  ionospheric  corn!  i  t  ions- - I  n  tin  previous  wink, 
it  was  assum'd  that,  for  daytime  c«>nd  i  t  i<*ns ,  l  h<  ilt  1 1  uii 


was  below  the  bottom  of  the  F-region  and  that  the  conduct ivit 
profile  in  the  neighborhood  of  h^  could  be  approximated  I"-  an 
exponential.  As  discussed  in  Section  III,  this  is  the  limit 
ht>  ■  ■  h^,  s  •  •  1  of  the  more  general  conductivity  profile.  It. 
this  limit, 


i(  2 


Since  ,  .  _  ,  1,  we  may  use  the  asymptotic  expansion  of  the 

•  function  to  obtain 


■  '  1 

.  ( 1  +  .  )  Ss  lo  j  - - -  -  log  2  -  log  s.  (  ’•2) 

Combined  with  by.  (!•*),  this  gives 

C’(  .  )  ■  hj  -  ;  l  l  (  -  )  -  2{  .  -  log  2)1  (  :■  .  ) 

where  ,  is  Ku  let's  constant.  For  the  i sot  topic  daytime 
ionospher-  of  Paper  1  ,  =  7  ,  so  that 

f  (  +  ,  _)  hj  +  '  2|  Y  +  2  (  ,  -  log  2 ) j  .  (  I  ) 

In  the  limit  s  1,  ,  ~  ,.  Thus,  apart  from  a  small  cotter 

t  ion  term,  which  wa.s  neglected  in  Paper  I,  t  he  e  i  tenvul  u<  • 
obtained  from  Fys.  (2'1)  and  (  ‘  4 )  ,  agrees  with  by.  (I). 

For  the  anisotropic  daytime  ionosphete  ot  Pa; it  11, 

Thus  , 
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and 
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( 3  S) 


whore  the  small  correction  term  has  been  neglected.  Since 
h  »  this  gives  the  same  eigenvalue  as  Eq .  (4)  . 

b.  Nighttime  ionospheric  conditions--It  was  assumed  in 
Paper  II  that,  for  nighttime  conditions,  a  sharp  E-rcgion 
bottom  was  encountered  before  the  altitude  h-^  was  attained, 
and  that  the  conductivity  increased  slowly  above  the  E-region 
bottom.  As  discussed  in  Section  III,  this  is  the  limit 
h  ■  h^,  s  >>  1  of  the  more  general  profile.  To  examine  this 

limit,  it  is  convenient  to  write 
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which  follows  from  Eqs.  (7),  (9),  and  (10)  with 
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In  the  limit  s  •  1,  |  |  1  and  .  (1  +  .•)  ~  ,  (1).  Furthermore, 


log  s  =  (h1  -  hF)  /  • 


(  3  K ) 


It  then  follows  from  Eqs.  (ri)  and  (  10)  that 
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To  lowest  order  in  e,  the  eigenvalue  obtained  from  Eqs. 
and  (39)  agrees  with  Eq.  (5). 


(40) 


(29) 


We  have  thus  demonstrated  that  the  results  of  Papers  I 
and  II  can  be  recovered  from  the  more  general  ionospheric 
conductivity  profile. 
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VI. 


APPROXIMATE  EXPRESSIONS  FOR  EXCITATION  FACTOR 


The  ELF  signal  strength  at  a  receiver  depends  not  only  on 
the  propagation  loss  between  transmitter  and  receiver,  but 
also  on  the  excitation  factor.  Predictive  calculations  of 
the  effects  of  ionospheric  disturbances  on  ELF  signal 
strengths  must  therefore  take  this  factor  into  account.  In 
Papers  I  and  II,  approximate  expressions  were  developed  only 
for  the  phase  speed  and  attenuation  rate.  In  this  section, 
the  theory  will  be  augmented  with  approximate  analytic  expres¬ 
sions  for  the  excitation  factor  as  well. 

For  the  calculation  of  the  excitation  factor,  it  is 
convenient  to  reformulate  the  theory  somewhat.  Following  a 
customary  procedure  in  ELF  calculations,  we  first  make  a 
Fourier  decomposition  in  the  transverse  plane  by  writing 


i3  ( i  z ) 


E  (  a  ,  v  ,  Z  ) 
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where  the  hat  is  used  to  denote  the  Fourier  transform.  We 
next  introduce  a  rectangular  set  of  ^-dependent  unit  vectors 
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and  resolve  all  Fourier  transforms  into  their  S-dependent 
components.  When  this  is  carried  out,  Maxwell's  equations 
separate,  in  the  isotropic  lower  ionosphere,  into  two  sets  of 
uncoupled  equations,  one  governing  the  propagation  of  TL  fields 
and  the  other  of  TM  fields.  The  equations  for  the  TM  fielus, 
which  include  the  TPM  mode,  arc 
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where  ?.  =  k  z  and 
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If  we  introduce  as  n  variable  the  wave  impedance 
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(4a) 


tne  TM  field  equations  can  be  combined  into  the  single  equation 
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winch  is  <i  nunlmmi,  first-order  differential  equation  el  tin 
Hiccnt  i  l  y pt  . 


2  0 


We  consider  a  small  horizontal  electric  dipole  of  strcnyt n 
IL,  oriented  along  the  x-direction  (;=0),  and  located  at  a 
height  zg  in  a  thin  uniform  "source  layer”  above  the  ground. 
The  solution  in  tne  source  layer  can  be  written  as 
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where  z  is  the  greater  of  z  and  z  and  z  is  the  smaller  of 
>  s 

these  two  quantities,  and 
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The  quantities  R  and  r  are  reflection  coefficients  to  be 
determined  from  the  boundary  conditions  at  the  upper  and  lower 
boundaries  of  the  source  layer,  respectively.  These  quant i t  m 
can  be  expressed  in  terms  of  the  wave  impedances  at  these 
boundaries  by  the  relations 
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where  we  have  placed  the  source  layer  at  tile  <u  c.und  <•'  0)  ana 

allowed  it  to  become  inf ini Lesima 1 ly  thin. 
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'idle  b‘  •unuary  conditions  require'  the  cont  mui  t  y  •. 
tandentiai  electric;  and  maanetic  field:;  at  all  inter: 


This  allows  us  to  evaluate  Z ^ ( 0 )  from  the  solutions  in  the 
yrounu.  In  deriving  the  approximate  eiyenvalues,  it  was 
assumed  that  the  ground  was  perfectly  conducting .  The  finite- 
conductivity  of  the  grounu  introduces  only  a  small  correction 
to  these  expressions.  however,  the  assumption  of  infinite 
ground  conductivity  cannot  be  made  in  deriving  the  excitation 
factor  for  a  horizontal  electric  dipole  at  the  ground,  since 
there  woulu  be  no  excitation  of  the  waveguide  at  all  under 
those  conditions.  The  solutions  in  the  ground  obeying  the 
radiation  condition  have  tire  form 
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where 
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It  then  follows  from  Eqs.  (4  ib)  and  (41)  that 
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The  impedance  Z^ (0)  is  determined  from  the  solutions 
above  the  source.  The  wave  equation  for*  the  impedance  can  b< 
written  as  the  integral  equation 
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the  contribution 


It  is  not  difficult  to  show  that,  for  7  ■  ■  h.  , 
to  the  integral  from  the  last  term  in  t.ne  square  brackets  is 
small.  ?nus , 


Z  =  ^(0)  +  z 


«  n  j ) 


where  nZ  is  given  by  Eg.  (44  j  .  Furthermore,  if  ,  little 

error  is  made  in  extending  the  range  of  integration  to  infin¬ 
ity.  The  wave  impedance  above  the  source  thus  takes  rile  for:-.. 


Z  -  V0)  +  V'  +  ko^o 

(h  '  x  '  h  1 )  (  .  4  ) 

where  ii  is  defined  by  the  equation 
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For  an  exponential  conductivity 
right-hand  side  of  Eg.  (28). 

From  tile  solutions  given  in 
near  h^,  the  impedance  somewhat 


profile,  hy  is  given  by  the 

Papers  I  and  II  for  the  regioi 
below  h  ^  takes  the  form 
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wher  e 


h^  =  £(l-  +  ,:_)  (  :J  7  ) 

and  £(.,:)  is  given  by  Eq.  (  50)  •  For  the  particular  profile 
treated  in  previous  work,  f  (.-,.•_ )  reduces  to  Eqs .  ( .>4)  ,  {  >:.)  , 

or  ( ) ,  as  the  case  may  be. 

By  requiring  the  leading  terms  of  Eqs.  (  r>4)  and  ( to 
agree  in  the  region  h  ■  where  both  are  valid,  we 

obtain 
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Z  ( 0 )  =  kQ(S  hQ  -  hx)  .  i  '  ■  ) 

a 

It  should  be  noted  that  at  this  point ,  we  may  regard  S“  as  an 
undetermined  parameter,  rather  than  the  TEM  eigenvalue. 

With  Z.  (0)  and  Z..(0)  now  determined,  we  can  express  the 
jj  U 

reflection  coefficients  r  and  R  in  terms  of  the  parameter  S  by 
means  of  Eqs.  (4'Ja)  and  (49b).  This,  finally,  allows  us  to 
express  E^  and  B^  in  terms  of  this  parameter.  The  expression 
for  B^  at  the  yrounu  becomes 
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where  wc  have  neglected  terms  of  relative  oi dor  n  ,  which 
are  small  at  ELI'. 

It  now  remains  to  carry  out  the  indicated  integration  in 
Eq.  (41a).  Carrying  out  the  integration  ovei  angle  in  S- space 
we  obtain 
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In  the  remaining  integration,  the  integral  can  be  extended  to 
-■  by  replacing  the  Bessel  function  by  an  appropriate  Hunkel 
function.  The  integral  can  then  be  evaluated  by  contour 
integration  in  the  complex  S-plune.  At  large  distances  from 
the  source,  the  primary  contribution  is  from  the  pole  of  the 
integrand  at 
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( o ; ) 


which  will  be  recognized  as  the  approximate  TL'M  eigenvalue 
obtained  previously.  This  may  be  considered  as  an  alternative 
derivation  of  that  expression.  From  the  residue  at  this  pole, 
we  obtain 
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whero  •  j  is  the  wavelength  in  live  space  and  we  ii.r.’c  useu 
large  argument  expansion  of  the  ilankel  function. 

We  del  me  tin  magnetic  excitat  ion  factor  ’  iry  tin 
equation 


i 

o 


2n  .  t 


^ . 


'O'  t) 


(<•  •:  ) 


which  differs  somewhat  from  other  definitions  (e.q.,  Bnnr.istc  r 
[  1975])  in  separatinq  the  effect  of  the  ionosphere  from  that 
of  the  ground.  With  this  definition,  th<  excitation  factor 
becomes 


,  _  -—l  ^-1/2  _  -3/4  —1/4 
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where  h.  is  given  by  Eq.  ( 17) .  Finally,  since  SqB  ,  the 

analogous  excitation  factor  for  the  vertical  electric  field  is 
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It  was  pointed  out  in  Paper  I  i  tiiat  Liu  m  r  i  ;•<  :.?  .i !  r  ac 
of  energy  flow  is  essentially  constant  up  to  an  altitude  ii  , 
above  which  it  falls  off  very  rapidly  with  altitude.  It  is 
tnerefore  not  surprising  that  the  effective  waveguide  height 
for  excitation  is  roughly  h  ,  rather  than  the  higher  reflect¬ 
ing  height  ii^. 


VII  . 


CONCLUSIONS 


The  approximate  theory  developed  in  Papers  I  and  II  for 
the  determination  of  ELF  propagation  constants  in  the  earth- 
ionosphere  waveguide  lias  been  extended  to  a  rather  more  gen¬ 
eral  ionospheric-conuuctivity  profile  than  those  considered 
previously.  The  results  of  Papers  1  ana  II  are  recovered  as 
limiting  cases  of  the  more  general  profile.  Simple  analytic 
expressions  have  also  been  derived  for  the  excitation  factor, 
which  is  another  important  parameter  in  ELF  propagation. 
Althougn  approximate  excitation  factors  have  not  been  compel eu 
with  full-wave  calculations,  it  is  reasonaule  to  assume,  on  thv 
basis  of  the  generally  excellent  agreement  between  approximate 
and  full-wave  propagation  constants,  that  they  represent  very 
good  approximations.  The  simple  analytic  expressions  derived 
in  this  work  allow  the  rapid  computation  of  pause  speeds , 
attenuation  rates,  and  excitation  factors  for  a  wide  range  of 
ionospheric  conditions  without  the  necessity  of  lengthy  full- 
wave  numerical  calculations .  The  theory  can  therefore  be 
quite  useful  in  the  study  of  the  effects  of  a  variety  of 
ionospheric  disturbances,  both  natural  and  artificial,  on 
LLF  communication  systems. 
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